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Abstract 

Recently, it has been shown, contrary to previous beliefs, that the k± 
distribution of quarks in a transversely polarized proton can be asymmetric. 
This "Sivers effect" had already been used to explain transverse single spin 
asymmetries (SSA) observed in inclusive pion production, p^ p — > -k X and 
p^ p — > ttX. In such channels, however, other mechanisms, like the "Collins 
effect" (a k± asymmetric fragmentation of a transversely polarized quark into 
pions), may generate SSA. The Sivers asymmetry is used here to compute 
SSA in Drell-Yan processes; in this case, by considering the differential cross- 
section in the lepton-pair invariant mass, rapidity and transverse momentum, 
other mechanisms which may originate SSA cannot contribute. Estimates for 
RHIC experiments are given. 



1. Introduction 



Single spin asymmetries in high energy inclusive processes are a unique testing 
ground for QCD; they cannot originate from the simple spin pQCD dynamics - dom- 
inated by helicity conservation - but need some non perturbative chiral-symmetry 
breaking in the large distance physics. Experiments are rich of puzzling single spin 
results and more are expected to come; new theoretical ideas and models have just 
appeared in the literature and a good understanding seems to emerge. 

Among the best known transverse single spin asymmetries (SSA) let us mention: 

• the old problem of the large polarization of A's and other hyperons produced 
in the scattering of unpolarized nucleons, pN — > A 1 * X [1]; 



• the large asymmetry 

da^ — da^ 



(1) 



observed for pion inclusive production, in p 1 * p — > n X and p^ p — > n X processes 

[2]; 

• the similar azimuthal asymmetry observed in semi-inclusive DIS processes, 
ttf -> InX [3]. 

All of these cannot be related to the elementary dynamics, but rather to non 
perturbative aspects of the nucleon and hadron structures. 

Among the several attempts [4] to explain the data within QCD dynamics we 
focus here on a phenomenological approach based on the generalization of the fac- 
torization theorem with the inclusion of parton intrinsic motion k± inside a nucleon 
and of hadrons relatively to the fragmenting parton. The cross-section for a generic 
process A B — ■> C X then reads: 

d(r=J2 fa/A(x a , k ±a )(g)f b/B (x b , k ±b )(g)da ab ^ c -(x a ,x b , k ±a , k ±b )(g)D c /c(z, k ±c ) (2) 

a,b,c 

where the /'s are the k± dependent parton distributions and the Z)'s the k± depen- 
dent fragmentation functions. 

The above QCD factorization theorem - with unintegrated k± dependent dis- 
tribution and fragmentation functions - has never been formally proven in general 
[5], but only for the Drell-Yan process (which is the issue of this paper) and for the 
two-particle inclusive cross-section in e + e~ annihilation (somehow a time-reversed 
Drell-Yan process) [6]. However, Eq. (2) has been widely used in the literature, 
starting from the pioneering work of Feynman, Field and Fox [7]. Actually, sev- 
eral papers have recently shown that intrinsic k±_'s are indeed necessary in order 
to be able to explain, within pQCD and the factorization scheme, data on large pr 
production of pions and photons [8]; without them the theoretical (collinear) com- 
putations would give results in some cases much smaller (up to a factor 100) than 
experiments. 
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When dealing with polarized processes the introduction of k± dependences opens 
up the way to many possible spin effects; these can be summarized, at leading twist, 
by new polarized distribution functions, 



A N 



fqtf = /,/pt (x, k±) - f q/p i (x, k±) = f q/p i (x, k±) - /,/pt (x, -k±) (3) 
fqt/p = f q y P (x, k ± ) - f q i /p (x, k±) = f q y p (x, k±) - f q y p (x, -k ± ) (4) 



and new polarized fragmentation functions, 



A^-Dft/qT = D h/q i (z, k ± ) - D h/q i (z, k ± ) = D h/q t (z, k ± ) - D h/q , (z, -k ± ) (5) 
& N D h y q = D h y q (z, fe ± ) - D h y q (z, k ± ) = D h y q (z, k ± ) - D h y q (z, -k ± ) (6) 



which have a clear meaning if one pays attention to the arrows denoting the polarized 
particles. All the above functions vanish when k± = and are naively T-odd. The 
ones in Eqs. (4) and (5), when written in the helicity basis, relate quarks of different 
helicities and are chiral-odd, while the other two are chiral-even. 

Similar functions have been introduced in the literature with different notations: 
in particular there is a direct correspondence [9] between the above functions and 
the ones denoted, respectively, by: [10], [11], and D^ T [10, 12]. 

In the recent comprehensive review paper on transverse quark polarization [13] 
the following notations are used: A N f q / p t = Aq/, A N f q y p = A°/ and A N D h / q i = 



The fragmentation in Eq. (5) is the Collins function [5], while the distribution 
in Eq. (3) was first introduced by Sivers [14]. 

Some of the above functions have been used for a phenomenological description 
of the observed SSA. Both Sivers [15] and Collins [16, 17] functions can explain the 
E704 data on An [2]: the Collins asymmetry might have been indirectly observed in 
HERMES data [3] and the polarizing fragmentation functions (6) can describe [18] 
the A polarization data [1]. 

Despite its successful phenomenology, the Sivers function was always a matter of 
discussions and its very existence rather controversial; in fact in Ref. [5] a proof of 
its vanishing was given, based on time-reversal invariance. Ways out based on initial 
state interactions [15] or non standard time-reversal properties [19] were discussed, 
but the situation remained uncertain and further phenomenology was performed 
making use of the Collins function only [16]. A similar criticism applied to the 
function in Eq. (4). 

Very recently a series of papers have resurrected Sivers asymmetry in its full 
rights. First, in Ref. [20], a quark-diquark model calculation including final state 
effects has given an explanation of the HERMES azimuthal asymmetry different 
from the one originating from Collins function. As a consequence, Collins recognised 
[21] that such a new mechanism is compatible with factorization and is due to the 
Sivers asymmetry (3); re-examining his original proof of the vanishing of A N f q / p i he 
finds it to be invalidated by the path-ordered exponential of the gluon field in the 
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operator definition of parton densities. The authors of Ref. [22] have confirmed the 
model results of [20]; a second paper by Brodsky, Hwang and Schmidt [23], along 
the same lines of the first one, shows that initial state interactions give rise to SSA 
in Drell-Yan processes. 

Some issues concerning factorizability and universality of these effects are still 
open to debate; however, we feel now confident to use Sivers effects - and equally 
all functions in Eqs. (3)-(6) - in SSA phenomenology. The natural process to test 
the Sivers asymmetry is Drell-Yan: in such a case there cannot be any effect in 
fragmentation processes and, by suitably integrating over some final configurations, 
also possible effects from transversely polarized quarks in an unpolarized proton, 
Eq. (4), do not contribute. SSA in Drell-Yan processes are particularly important 
now, as ongoing or imminent experiments at RHIC will be able to measure them. 

The paper is organized as follows: in Section 2 we present the formalism necessary 
to discuss single spin asymmetries in Drell-Yan processes, taking into account the 
transverse motion of partons in nucleons. Some simple analytical formulae are also 
derived, assuming simplified functional k± dependences. In Section 3 we give some 
numerical estimates for SSA at RHIC, while comments and conclusions are gathered 
in Section 4. 



2. SSA in Drell-Yan processes, formalism 

Let us consider a Drell-Yan process, that is the production of pairs in the 
collision of two hadrons A and B; there is no need of any fragmentation function 
and Eq. (2) reads: 



da = J2 J [ dx a d 2 k ±a dx b d 2 k ±b f a /A(x a , k ±a ) f h / B {xb, k±b) da ab ~ 

ab 



(7) 



and, from the Sivers asymmetry of Eq. (3), the difference between the single trans- 
verse spin dependent cross-sections da^ for B — > i + 1~ X and da 1 for A^ B — > 
£+ £~ X is 



ab^e+r 



da^ - da 1 = Y^ J \dx a d 2 k ±a dx b d 2 k ±b A N f a/A i (x a , k ±a ) f b / B (x b , k± b ) da 

ab 



The elementary cross-section da for the process a{p a )b{p b ) — > £ + (p + )£ is 
given by: 



1 d 3 p + d 3 p. 
da = — 



S 4 (p a + Pb ~ P+ ~ P-) | M ab ^ e+e - 



2s 2E + 2£„ (2tt) 2 
We consider the differential cross-section in the variables 

1 , qo + QL 



s = M 2 = (p a +p b y = q 2 



y 



In 

2 g - qL 



(9) 



(10) 
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that is the squared invariant mass, the rapidity and the transverse momentum of the 
lepton pair; q , q T and qi are respectively the energy, transverse and longitudinal 
components, in the A-B center of mass frame, of the four- vector q = p a +Pb = P++P-- 
Using the relations: 

f^ = d A p_5(p 2 _) P- = q~P + dM 2 dy = 2dq dq L , (11) 

Eq. (8) can be written as 

dV dV 1 



^2 / [dx a d 2 k ±a dx b d 2 k± b 



dydM 2 d 2 q T dydM 2 d 2 q T 2 c , 
A N f a/A r{x a , k ±a ) f b/B {x b , k ±b ) 5\p a +p b - q) of (12) 

where af is the total cross-section for the ab — > £ + l~ process: 

d 3 p+ 1 1 
2E^ (2tt) 2 2M 2 

Analogously 



< = / l£ (2W '«« - P +> 2 > I M — «> I' ' (13 > 



dV dV 

dydM 2 d 2 q T + 'd~ydM 2 d 2 q T 



= Yl J [ dXa d2 k^a dx b d 2 k ±b 



fa/A(x a ,k ±a ) f b / B (x b ,k ±b )5 (Pa+Pb-q)^ ( 14 ) 

which is twice the unpolarized cross-section. 

For Drell-Yan processes the dominating electromagnetic elementary interaction 
is qq — > 7* — > so that a,b = q,q with q — u,u, d, d, s, s and: 

Eqs. (12) and (14) allow to compute a single spin asymmetry A^v, Eq. (1), 
for Drell-Yan processes and for differential cross-sections measuring the lepton pair 
invariant mass M, rapidity y and transverse momentum q T ; notice that we do 
not look at the angular distribution of the lepton pair production plane, which is 
integrated over. 

Let us now fix in more details our kinematical configuration. We take the hadron 
A as moving along the positive z-axis, in the A-B cm. frame and measure the 
transverse polarization of hadron A, Pa, along the y-axis, as shown in Fig. 1. 
Neglecting masses, the four-momenta of hadrons and partons are: 

^4 = ^(1,0,0,1) p B = ^ (1,0 ,0,-1) (16) 



Pa = * fl ^ll + %%, I" V) (17) 
2 \ X 2 S X a ^S X 2 S J 
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f k\ h 2k Lb 
Pb = x b — 1 + — , r , 

2 \ x b s x^s 



-1 + ^ 

xis , 



(18) 



with 



q = Pa + Pb = (go, <?t> 9l) = (yM 2 + g| coshy, q T , ^ M 2 + q 2 , sinhy) , (19) 

where the lepton pair rapidity y and invariant mass M are defined in Eq. (10). 

The four-momentum conservation S function of Eqs. (12) and (14) contains the 
factors: 



5{E a + E b - g ) S(p za + p zb - q L ) = 



s((x a - Xb)^Y 



k 2 h 2 

x a s x b s 



k 2 



k 2 
^Ib 



x a s x b s 

We shall consider kinematical regions such that: 

q 2 T « M 2 k\ a b ~ q 2 T , 

where Eq. (20) simplifies into the usual collinear condition: 



- 5{E a + E b - q )S{p za + p zb - q L ) 



\ A 



M 



M 



(20) 



(21) 



(22) 



In such a region [q 2 ^ <C M 2 , qx — k±] then one has 

E 9 e 2 J d 2 k ±q d 2 k ± g 5 2 (k ±q + k ±q - - q T ) A N f q/Al (x q , k ±q ) f q / B (x q , fe±g) 
2E 9 e2J 'd 2 k ±q d 2 k ±q S 2 (k ±q + k ±q - q T ) f q / A (x q , k ±q ) fq/ B (x q , k± q ) 



A 



N 



with x q and x q fixed by Eq. (22) with a,b = q,q and q = u,u, d, d, s, s. 



(23) 



2.1 Other mechanisms for SSA in Drell-Yan processes 

Before discussing further Eq. (23), let us comment on other possible origins of 
SSA. Let us consider first the SSA generated by the distribution function in Eq. (4), 
as compared with Sivers mechanism, Eq. (3), which we are considering here. The 
main point is that, as we have seen, Eq. (3) gives a contribution to A N of the type: 



AN fq/Al (Xa, k± a ) ® fq/B{x b , k ±b ) <g> do™' 



(24) 



where da is the elementary unpolarized cross-section, while Eq. (4) leads to a con- 
tribution of the kind [11] 



E M^a, *!±a) ® A N f q y B ( Xb , k ±b ) <g> dA&^ 



(25) 
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where h\ q is the transversity of quark q (inside hadron A) and dAa is the double 
transverse spin asymmetry da^ — da^. Such an elementary asymmetry has a cos 20 
dependence [11, 24], where is the angle between the transverse polarization di- 
rection and the normal to the £ + £~ plane; when integrating over all final angular 
distributions of the £ + £~ pair - as we do by looking only at variables (10) - the 
contribution of Eq. (25) vanishes. 

There exist in the literature other mechanisms to generate SSA in Drell-Yan 
processes [25, 26, 27, 28], based on higher twist quark-gluon correlation functions in 
a generalized pQCD factorization theorem [29]. All of them lead to expressions of 
An depending on the angles between the polarization direction and the final lepton 
pair plane [30], which require observation of these angles to be detected and vanish 
upon integration. 

In the model of Ref. [31] a SSA for the differential cross-section (14) is computed: 
in that model a non vanishing asymmetry is achieved by introducing orbital angular 
momentum and surface effects in the distribution of valence quarks, resulting in 
somewhat ad hoc correlations between their polarization and k± distribution. It 
might be that the mechanism is somehow related to the Sivers effect, although the 
issue should be further investigated. 

In Refs. [20] and [23] - which are meant to be and provide important pedagogical 
examples of SSA - the processes considered are somehow academic and related by 
crossing. They are respectively 7*p^ — > q(qq)o and qp^ — > 7* (qq)o where (qq)o is 
a spectator scalar diquark. The SSA obtained in these two cases turn out to be 
opposite, as predicted by Collins [21]. 

2.2 A simple analytical model 

In order to give numerical estimates in the next Section, we introduce here a 
simple model for the Sivers asymmetry (3), and for the unpolarized distributions, 
which is similar to the one introduced for the polarizing fragmentation function in 
Ref. [32] and has the advantage of giving analytically integrable expressions for A N . 
Such a model shows explicitely how the asymmetry originates and depends on M, 
y and q T . 

Let us start from the most general expression for the number density of quarks 
q, inside a proton with transverse polarization P and three-momentum p; the quark 
has a transverse momentum k± and its polarization is not observed. One has 

fq/ P i(x, k±) = f q /p(x, k±) + ^ A N f q/p1 (x, k±) P ■ p x k ± , (26) 

consistently with Eq. (3) 

/,/pt (x, k ± ) + f q/pl (x, k ± ) = 2 f q/p (x, k ± ) = 2 f q/p (x, k ± ) (27) 
f q /pi (x, k±) - f q/p i (x, k±) = A N f q/ ,^ (x, k±) = A N f q/p i (x, k±) P • p x k ± . (28) 

With the configuration of Fig. 1 one simply has P ■ p x kj_ = (kj_) x = cos0fc ± . 
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The Sivers function A N f q / p i(x, k±) must obey the positivity bound: 

l AiV /g/pT (x, kj 



< 1 



2fq/ P (x, k±) 

We consider simple factorized and Gaussian forms: 

(3 2 2 2 

f q/p (x, k ± ) = f q / p (x) g(k ± ) = f q / p (x) — e - ^ fe ' , 



(29) 



A JV / 9 / P T(x,A; ± )=A w / g / P T(a;)MA; ± ). 



(30) 
(31) 



Notice that j3 can depend on x; we have assumed the same kj_ dependence for all 
flavours, as intrinsic k± of quarks originates via confinement and gluon emissions, 
which should be flavour independent processes. Notice also that [3 2 = 1/ (k±) and 
that fd 2 k ± f q / p (x, kj_) = fq/ p (x). 

In order to obviously satisfy the bound (29) we write 



A N f q/p ,(x)=2Af q (x)f q/p (x) 

h(k ± )=H(k ± )g(k ± ). 
so that Eq. (29) simply becomes 

\Nq(x)H(k ± )\<l Vx,k ± 



(32) 
(33) 

(34) 



We actually impose Af and H to be separately smaller than 1 in magnitude, 
by choosing simple analytical functional forms and dividing each of them by its 
maximum value: 



Mq(x) 



N n x a "(l-x) 



a? bq" 



\Nq\ < 1 



-(a 2 -I3 2 )ki , a> P 



7~t(kj_) = ^2e{a 2 -p 2 ) k± exp 
Eqs. (30)-(33) and (36) imply: 

A iY / 9/p T(x, k ± ) = 2Af q (x) f q/p (x) F ^J2e(a 2 -(3 2 )k ± e- a2 



(35) 
(36) 

(37) 



Inserting the above choice of A N f(x, k±) and f(x, k±) into Eq. (23) one can 
perform analytical integrations. The numerator of Eq. (23) results as: 



N(A N ) 



3/54 



7T (a 2 + p 2 y 



2e 



a 2 - p 

~T 2 



2\ V2 



(q T ) x exp 



a 2 (3 2 2 
Qt 



a 2 + p 2 



(38) 
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where a = a(x q ), f3 = /3(x q ), f3 = f3(xq) with x q = Me 9 /VS) x q = M e y '/y/s. The 
denominator is 



D(A N ) 



1 (3 2 P 2 
Ti (3 2 + p- 



i ex P 



I3 2 p 2 



^QT 2 Y,<%fq/p(Xq)fq/p(x 9 ), (39) 



where again f3 stands for f3(x q ) and f3 for (3(x q ). 
The asymmetry (23) in this simple case reads: 




) 



1 E q e 2 q A N f q/p r(x q )f q/p (x q ) _ 

2 J2 q 6 q fq/p( X q) fq/p( X q) 



(40) 



For a generic configuration, different from the one of Fig. 1, one simply replaces 



When, as it is often the case in the literature, a and f3 are taken to be independent 
of x, Eq. (40) simplifies to: 



The same simplified expression holds also at y — 0, where one has x q = x q = 



3. SSA in Drell-Yan processes, numerical estimates 

In this Section we shall present numerical estimates of the SSA for the Drell- 
Yan process, originating from Sivers effect, in kinematical configurations relevant 
for RHIC experiments. 

We follow the simple analytical model discussed in the previous Section; we need 
to fix the parameters of the model, that is the functions /3(x), a(x) and the flavour- 
dependent coefficients N q , a q , b q appearing respectively in Eqs. (30), (37) and (35). 
The unpolarized and k± integrated partonic distributions are chosen among the sets 
available in the literature; in particular, we adopt the GRV94 set [33]. 

Some information on /3, entering the dependence in the unpolarized parton 
distributions, can be obtained by considering the available data on the unpolarized 



(q T ) x with P ■ p x q T . 




(41) 
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cross-sections for inclusive particle production processes. This not only supplies 
information on f3, but allows a crucial consistency check of our formalism, as it 
shows whether or not we can reproduce the unpolarized data. 

Experimental data are presently available for pion, prompt photon production 
and for Drell-Yan processes, in different kinematical configurations. As mentioned 
in the introduction, several papers [8] have studied in details the unpolarized cross- 
sections for these processes within perturbative QCD at LO and NLO, emphasizing 
the crucial role of intrinsic transverse momentum. 

We have performed an independent attempt of reproducing the available data 
within the same approach utilized in the calculation of the SSA. Since SSA have 
been studied at present only at leading order and leading twist, the same level of 
accuracy has been adopted for the unpolarized cross-sections, keeping in mind that 
there might be NLO corrections (the so-called K factors) the value of which may 
vary, depending on the process and on the kinematical configuration considered, 
approximately between a factor 1 and 3. 

Our aim is not to reproduce as accurately as possible the unpolarized cross- 
sections, as in [8], but rather to show that within the same approach used for the 
SSA it is possible to reproduce the unpolarized cross-section values, up to an overall 
factor between 1 and 3, due to NLO corrections, scale-dependences, etc. These 
corrections should have little effect on the SSA, largely canceling in the ratio of 
cross-sections. 

A full account of the combined study of unpolarized cross-sections and SSA in 
different processes is outside the scope of this paper, and will be given in a separate 
forthcoming paper [34]. It suffices to state here that, as a result of this comprehensive 
analysis, a value of (3 independent of x, (3 — 1.25 (GeV/c) -1 , corresponding to 

\J (kj_) = 0.8GeV/c, allows a good description of the unpolarized processes and is 
appropriate for our scope. This value is in good agreement with the results of the 
abovementioned papers devoted to unpolarized cross-section calculations [8], and 
will be adopted in the sequel. 

One might argue that the value of (3 varies with energy and that RHIC works at 
cm. energies much larger than most of previous experiments; this could be tested 
when RHIC data on unpolarized cross-sections will become available and, for the 
moment, we keep using the value (3 = 1.25 (GeV/c) -1 , which corresponds to an 
already large ( k± ) value. We only notice that in our approach changes in (3 should 
influence only the dependence of the SSA on |qf T |. In the following, we will show 
such a dependence by presenting some of our results for (3 = 0.83 (GeV/c) -1 , which 

means yf{ k\ ) = 1.2 GeV/c. 

Adopting a f3 parameter independent of x and, as we do, of quark flavour, sim- 
plifies our analytical expressions as one gets f3 = (3, which allows the use of Eq. (41). 

As a second step, we must fix the parameter a determining the gaussian shape 
of the Sivers function, see Eq. (37). Since the positivity bound, Eq. (29), requires 
a > /3, one can write a 2 = fi 2 /r, where r is a numerical parameter in the range (0, 1). 
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A similar procedure was already followed by us in a previous work on transverse 
hyperon polarization [32], and an optimal guess for the parameter r was found, 
r ~ 0.7, value which we also adopt here. The same value allows a good description, 
using Sivers function, of the E704 data on [2, 34]. 

The next step is to obtain information on the x-dependent part of the Sivers 
function, A N f q / p t(x). Indications on the values of the parameters entering Eq. (35) 
can be obtained by exploiting the data on SSA measured by the E704 experiment 
[2] for the process p^ p — > ir X. 

Estimates of the Sivers function from E704 results were already given and used to 
make predictions for SSA in similar processes in [15]. In those first papers, however, 
a number of simplifying approximations were adopted which are not kept in the 
present, more general, analysis. Those assumptions amounted to consider k± effects 
only minimally, in places where their neglect would lead to vanishing results; also, 
no gaussian k± distribution was introduced, but rather a simple <5-like dependence 
which leads to the effective use of an average k± in a simplified (planar) kinematics. 

Although the simplified procedure might be a good approximation for single spin 
asymmetries, which are ratios of cross-sections, a more general and refined analysis, 
using full kj_ dependences according to Eq. (2), is in progress and a detailed analysis 
will be presented elsewhere [34]; we only state here that good reproductions of both 
the cross-sections and the E704 asymmetry are possible and a good set of values for 
the parameters N q , a q and b q of the Sivers function is: 

N u = 0.5 a u = 2.0 b u = 0.3 

N d = -1.0 a, = 1.5 b d = 0.2 . (42) 

Notice that we are assuming that only valence quarks in the proton give a nonvan- 
ishing Sivers function. 

The above values allow a good description of the E704 data on the SSA in the 
p^ p — > n X process [2], assuming that A N is totally generated by Sivers asymmetric 
distribution. We use them here in order to estimate A N for a D-Y process at a 
much higher energy: we are well aware that this might be a rough extrapolation, as 
we neglect pQCD evolution, other possible mechanisms, and so on. However, since 
information on Sivers function is so scarce, even rough estimates of the SSA in Drell- 
Yan processes are important and useful, and we shall proceed with our program, 
keeping in mind all these remarks. 

One further uncertainty concerns the sign of the asymmetry: as noticed by 
Collins [21] and checked in Ref. [23], the Sivers asymmetry has opposite signs in 
Drell-Yan and SIDIS, respectively related to s-channel and t-channel elementary 
reactions. As in p — p interactions we expect, at large xp, a dominant contribution 
from t-channel quark processes, we think that the Sivers function extracted from 
p — p data should be opposite to that contributing to D-Y processes. Our numerical 
estimates will then be given with the same parameters as in Eq. (42), changing the 
signs of N u and Nj. Given these considerations, even a simple comparison of the 
sign of our estimates with data might be significant. 
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We have thus completed the choice of parameters of our model, and we are now 
able to give predictions for the Drell-Yan process, without forgetting all the cautious 
comments made above. 

We only add that a different modelization of the Sivers function, originally pro- 
posed by Collins for his function [5] and often adopted in the literature, could be 
used. This functional form does not allow to perform exact analytical integrations 
on transverse momentum, and we prefer using the form proposed here. However, we 
have checked that for consistent choices of the parameters, the two parameteriza- 
tions give very similar results, since they mainly differ at relatively large transverse 
momentum, where the gaussian dumping is already effective. 

An important point has to be noticed regarding the parameterization of the 
Sivers function. In order to reproduce the E704 results, it has to have a valence-like 
behaviour. That is, A N f q / p i(x)/ f q / p (x) ~ x a , where a > 0, for small x. 

This is a very general feature, somehow required by the experimental data, which 
may have strong effects for the SSA in Drell-Yan processes, leading to tiny values of 
A N when the x values explored are relatively small. This prediction and its extent 
of validity can be tested at RHIC. In what follows, we shall show some estimates 
for the SSA in realistic kinematical regions: 

6 GeV < M < 10 GeV; -2 < y < 2 . (43) 

To show the dependence on the range of invariant mass covered, we will consider 
also the alternative set of kinematical cuts 

10 GeV < M < 20 GeV; -2 < y < 2 . (44) 

These kinematical regions should be easily explored in RHIC experiments, while 
assuring that the Drell-Yan process is the dominating contribution to lepton pair 
production. Notice that with the above cuts the lowest possible value of x reached 
is around 10~ 3 . 

In our simple analytical model, the dependences of on q T and x q , x q (or M, 
y) are completely uncorrelated, and may be treated separately. In fact, assuming 
(3 = (3 and a 2 = (3 2 /r, Eq. (40) now reads 



A N (M,y,q T ) 



= 2- 



A(M,y) 



1 



r 



2e- 



1 



1/2 



P q T cos (f)q T exp 



1 1 



2 1 + r 



X 



1 Eg A N f g/p1 ( Xq ) fq/p(Xg) 

2 J2q e q fq/p( X q) fq/p( X q) 



(45) 



where <p q is the azimuthal angle of defined in Fig. 1. 

The second line of the above equation, the function Q(q T ), describes the q T 
dependence of the asymmetry, which turns out to be very simple and directly related 
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to the k± behaviour of the Sivers function (via the parameter r). One can easily 
see that Q(qr) has a maximum when 

q T = q™ = ^/(l + r )/(l-r)//? , (46) 

where its value is Q{qx) = Qm = [2r/(l + r) ] 3 / 2 . Notice that the position of the 
maximum depends on the parameter (3, while Qm only depends on r = f3 2 /a 2 . In 
particular, when r = 0.7, q^ ~ 2.38//3, and Q M ^ 0.75. 

In order to collect statistical significance, one could integrate the asymmetry over 
q T , but this would lead to a vanishing result, due to the cos0 9t factor. However, 
one can think of integrating over <p qT in the range [0, 7r/2] only, or alternatively of 
taking into account the change of sign in the different quadrants. In both ways, the 
4>q T integration gives an overall factor / d<p | cos0|/ J d<p = 2 /it. 

We can then consider the magnitude of the asymmetry averaged over q T up to 
an upper value of qx = qri- In our simple model, Eq. (45), one obtains 



r ^f 2e ^) 1/2 {V^^H-^^e- w2 } (i_ e -(i+0» 2 /2) 1 A(M,y) 



_ P T1 d 2 q T | A N (M, y, q T ) \ da _ J qT1 d 2 q T (da^ - do±) 
{ | A N {M,y, q T ) | ^ - jqri ^ ^ - jqTi ^ ^ + ^ 

ix \ 1 + r 

= Q(q T1 )A(M,y) (47) 

where da stands for da/dydM 2 d 2 q T , w = f3 ?ti/a/1 + r and A(M,y) is given by 
the last line of Eq. (45). In particular, when qx\ — > oo, one finds the very simple 
result: 

/2e\ 1 / 2 , /l-rX 1 / 2 
( I A N (M, y, q T ) I )oo = ( - ) ^' V2 ( — r ) A ^ V) ■ ( 48 ) 

This expression holds with a good accuracy already at qxi > 1//3. 

Notice again that, concerning the intrinsic motion dependence, the above average 
is independent of (3 and depends only on r, via a function which has a maximum at 
r = (vl~0 — l)/3 ~ 0.72, very close to the value r = 0.7 adopted in our analysis. 

In Figs. 2 to 5 we show our model estimates for the SSA, starting from Eq. (45). 
In Fig. 2 we plot as a function of y: the asymmetry is averaged over M, in the 
two kinematical ranges 6 < M < 10 GeV and 10 < M < 20 GeV, corresponding to 
the solid and the dashed curve respectively. We have fixed qr = q^f ■, Eq. (46), and 
<pq T = 0, which maximizes the q T -dependent part of the asymmetry. 

In Fig. 3 we show the same asymmetry, averaged over different rapidity ranges, 
\y\ < 2 (solid curve) and < y < 2 (dashed curve). 

We also show the dependence of An on xf in Fig. 4 and on x q in Fig. 5; these 
partonic variables ( q ) are obviously related to y and M via Eq. (22). The 

integration over the non observed variables is done in such a way as to correspond to 
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the kinematical ranges — 2 < y < 2, 6 < M < 10 GeV (solid curve) and — 2 < y < 2, 
10 < M < 20 GeV (dashed curve). 

All asymmetries in Figs. 2-5 are evaluated at qr = qp 1 : the value of the SSA at 
different qp values may be obtained by simply rescaling the results of Fig.s 2-5 by 
the factor Q{qp) / Q{qp ) ■ This factor is plotted in Fig. 6. In order to show how its 
behaviour depends on the parameter f3, we present results for (5 = 1.25 (GeV/c) -1 
and = 0.83(GeV/c)- 1 . 

Analogously, one can obtain the magnitude of the asymmetry averaged over q T , 
Eq. (47), by rescaling the values given in Figs. 2-5 by the appropriate factor: this 
factor is obtained by dividing Q(qri), see Eq. (47), by Q(qp), and it is shown in 
Fig. 7. 

Our numerical estimates show that A N can be well measurable within RHIC 
expected statistical accuracy. The actual values depend on the assumed functional 
form of the Sivers function and its role with valence quarks only. This reflects in 
the increase of to sizeable values at large x q , Fig. 5, and at large xp, Fig. 4; in 
fact, xp large and positive implies a large x q > xp. 

4. Comments and conclusions 

The single transverse spin phenomenology, within QCD dynamics and the fac- 
torization scheme, is a rich and interesting subject. It combines simple pQCD spin 
dynamics with new long distance properties of quark distribution and fragmenta- 
tion; the experimental measurements are relatively easy and clear, many have been 
and many more will be performed in the near future, both at nucleon-nucleon and 
lepton-nucleon facilities. 

Very recently a large single transverse spin asymmetry has been observed at 
RHIC, at the very first spin measurement at y/s = 200 GeV, in p T p — > tt X processes 
[35]; despite the large energy the asymmetry is not negligible, contrary to naive 
expectations. 

The approach adopted here - and in many previous papers - requires the explicit 
control of parton intrinsic motion, which cannot be simply integrated over in the 
different non perturbative functions and the elementary dynamics, but must be kept 
into account wherever it can give new effects. Actually, this is true also for a correct 
computation of the unpolarized cross-sections, as it has been known for a long time, 
although somewhat forgotten. When dealing with spin dependences, the intrinsic 
motion is a rich and unexpected source of many new effects. 

We have presented here the explicit formalism for computing single transverse 
spin asymmetries in Drell-Yan processes, within a generalized QCD factorization 
theorem formulated with dependent distribution functions. Simple gaussian 
forms have been assumed and available data from other processes have been ex- 
ploited, in order to give estimates for single spin effects in D-Y production at RHIC, 
which should be of interest for the incoming measurements. Again, sizeable and 
measurable values have been found. 
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Our approach can and will be extended to the study of any A B —> C X process, 
at large energy and moderate to large momentum transfer; the parton intrinsic 
motion is relevant both in polarized and unpolarized processes and the resulting 
phenomenology might explain or anticipate many subtle and unexpected results. 
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Figure captions 

Fig. 1: The kinematical configuration considered in this paper. The 7* four- 
momentum defines all our observables; the dependence on the angle between the 
7*-^ and the 7*-(£ + £~) planes is integrated over in Eq.s (12) and (14). 

Fig. 2: The single spin asymmetry An for the Drell-Yan process, see Eq. (45), at 
RHIC energies, y/s = 200 GeV, as a function of the rapidity y and averaged over the 
invariant mass ranges 6 < M < 10 GeV (solid line) and 10 < M < 20 GeV (dashed 
line). Results are given at qp = q? , see Eq. (46), and <pq T = 0, which maximizes 
the effect. Furthermore, we have used r = 0.7, the parameters of Eq. (42) for the 
Sivers function (see text for further details) and the parameterization GRV94 [33] 
for the unpolarized parton distributions. Notice that the asymmetry is practically 
negligible in the range y < 0. 

Fig. 3: The single spin asymmetry An for the Drell-Yan process, see Eq. (45), at 
RHIC energies, y/s = 200 GeV, as a function of the invariant mass M and averaged 
over the rapidity ranges — 2 < y < 2 (solid line) and < y < 2 (dashed line). 
Results are given at qp = qx ■> see Eq. (46), and <P Qt = 0, which maximizes the 
effect. Furthermore, we have used r = 0.7, the parameters of Eq. (42) for the Sivers 
function (see text for further details) and the parameterization GRV94 [33] for the 
unpolarized parton distributions. 

Fig. 4: The single spin asymmetry An for the Drell-Yan process, see Eq. (45), at 
RHIC energies, y/s = 200 GeV, as a function of the Feynman variable Xp — Xq — Xq 
and averaged over the rapidity and the invariant mass in the ranges —2<y<2, 
6 < M < 10 GeV (solid line) and -2 < y < 2, 10 < M < 20 GeV (dashed line). 
Results are given at qp = qp , see Eq. (46), and <pq T = 0, which maximizes the 
effect. Furthermore, we have used r = 0.7, the parameters of Eq. (42) for the Sivers 
function (see text for further details) and the parameterization GRV94 [33] for the 
unpolarized parton distributions. The two curves almost coincide but the solid line, 
corresponding to a lower invariant mass range, cannot reach values of xf <^ 0.36 
within the given y range. Notice that the asymmetry is practically negligible in the 
range xp < 0. 

Fig. 5: The single spin asymmetry A N for the Drell-Yan process, see Eq. (45), at 
RHIC energies, ^fs = 200 GeV, as a function of x q and averaged over the rapidity 
and the invariant mass in the ranges —2 < y <2, 6<M<10 GeV (solid line) and 
-2 < y < 2, 10 < M < 20 GeV (dashed line). Results are given at q T = qp 1 , see 
Eq. (46), and 9t = 0, which maximizes the effect. Furthermore, we have used r = 
0.7, the parameters of Eq. (42) for the Sivers function (see text for further details) 
and the parameterization GRV94 [33] for the unpolarized parton distributions. The 
two curves almost coincide but the solid line, corresponding to a lower invariant 
mass range, cannot reach values of x q > 0.37 within the given y range. 
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Fig. 6: The factor Q(q T )/Q(q¥)i see Eq.s (45) and (46), plotted function 
of q T , for = 1.25 (GeV/cT 1 (solid line) and (3 = 0.83 (GeV/cT 1 (dashed line), 
corresponding respectively to ( k\ ) 1 / 2 = 0.8 GeV/c and ( k\ ) x l 2 = 1.2 GeV/c. This 
factor can be used to rescale the asymmetries given in Fig.s 2-5, at qr = q^ \ to 
their values at qx different from q^f. 

Fig. 7: The factor Q(qTi)/ Q(Qt)i see Eq.s (47) and (46), plotted as a function 
of q T1 , for = 1.25 (GeV/cT 1 (solid line) and p = 0.83 (GeV/c)" 1 (dashed line), 
corresponding respectively to ( k\ ) 1 / 2 = 0.8 GeV/c and ( k\ ) 1 / 2 = 1.2 GeV/c. This 
factor can be used to obtain from the asymmetries given in Fig.s 2-5 (at fixed q T , 
q T = q^ an d 4>q T = 0) the corresponding asymmetries averaged over q T up to 
\Qt\ = Qti (see text for further details). 
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Figure 1: The kinematical configuration considered in this paper. The 7* four- 
momentum defines all our observables; the dependence on the angle between the ^*-z 
and the ^f*-{l + l~) planes is integrated over in Eq.s (12) and (14). 
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Figure 2: The single spin asymmetry An for the Drell-Yan process, see Eq. (45), at 
RHIC energies, y/s = 200 GeV, as a function of the rapidity y and averaged over the 
invariant mass ranges 6 < M < 10 GeV (solid line) and 10 < M < 20 GeV (dashed line). 
Results are given at qx = q? see Eq. (46), and <p q = 0, which maximizes the effect. 
Furthermore, we have used r = 0.7, the parameters of Eq. (42) for the Sivers function (see 
text for further details) and the parameterization GRV94 [33] for the unpolarized parton 
distributions. Notice that the asymmetry is practically negligible in the range y < 0. 
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Figure 3: The single spin asymmetry An for the Drell-Yan process, see Eq. (45), at 
RHIC energies, y/s = 200 GeV, as a function of the invariant mass M and averaged over 
the rapidity ranges — 2<y<2 (solid line) and < y < 2 (dashed line). Results are given 
at qr = q? \ see Eq. (46), and <p q = 0, which maximizes the effect. Furthermore, we 
have used r = 0.7, the parameters of Eq. (42) for the Sivers function (see text for further 
details) and the parameterization GRV94 [33] for the unpolarized parton distributions. 
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Figure 4: The single spin asymmetry An for the Drell-Yan process, see Eq. (45), at 
RHIC energies, y/s = 200 GeV, as a function of the Feynman variable and 
averaged over the rapidity and the invariant mass in the ranges — 2<y<2,6<M<10 
GeV (solid line) and -2 < y < 2, 10 < M < 20 GeV (dashed line). Results are given at 
It = Qt i see Ecf (46), and (j)g T = 0, which maximizes the effect. Furthermore, we have 
used r = 0.7, the parameters of Eq. (42) for the Sivers function (see text for further details) 
and the parameterization GRV94 [33] for the unpolarized parton distributions. The two 
curves almost coincide but the solid line, corresponding to a lower invariant mass range, 
cannot reach values of xf > 0.36 within the given y range. Notice that the asymmetry is 
practically negligible in the range xp < 0. 
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Figure 5: The single spin asymmetry An for the Drell-Yan process, see Eq. (45), at 
RHIC energies, y/s = 200 GeV, as a function of x q and averaged over the rapidity and the 
invariant mass in the ranges — 2 < y < 2, 6 < M < 10 GeV (solid line) and — 2 < y < 2, 
10 < M < 20 GeV (dashed line). Results are given at qr = Qt > see Eq- (46), and (j)q T = 0, 
which maximizes the effect. Furthermore, we have used r = 0.7, the parameters of Eq. (42) 
for the Sivers function (see text for further details) and the parameterization GRV94 [33] 
for the unpolarized parton distributions. The two curves almost coincide but the solid line, 
corresponding to a lower invariant mass range, cannot reach values of x q > 0.37 within 
the given y range. 
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Figure 6: The factor Q{qr) / Q{qr ) ■, see Ec l- S an d (46), plotted as a function of qT, 
for P = 1.25 (GeV/c) _1 (solid line) and ft = 0.83 (GeV/c) -1 (dashed line), corresponding 
respectively to {k\ ) l l 2 = 0.8 GeV/c and {k\ ) 1 / 2 = 1.2 GeV/c. This factor can be used 
to rescale the asymmetries given in Fig.s 2-5, at qr = q^ , to their values at qx different 
from qjf . 
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Figure 7: The factor Q(qri) / Q(Qt ') > see Eq-s (47) and (46), plotted as a function of qri, 
for (3 = 1.25 (GeV/c) -1 (solid line) and f3 = 0.83 (GeV/c) _1 (dashed line), corresponding 
respectively to {k\ ) 1/2 = 0.8 GeV/c and {k\ ) 1/2 = 1.2 GeV/c. This factor can be used 
to obtain from the asymmetries given in Fig.s 2-5 (at fixed q T , qx = q? and </>g T = 0) 
the corresponding asymmetries averaged over q T up to \q T \ = qxi (see text for further 
details). 
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